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Abstract. Two particle correlations are used to extract information about the characteristic 
size of the system for proton-proton collisions at y / s=900 GeV measured by the ALICE (A 
Large Ion Collider Experiment) detector at CERN. The correlation functions obtained show the 
expected Bose-Einstein effect for identical particles, but there are also long range correlations 
present that shift the baseline from the expected flat behavior. A possible source of these 
correlations is the conservation of energy and momentum, especially for small systems, where 
the energy available for particle production is limited. A new technique, first introduced by 
the STAR collaboration, of quantifying these long range correlations using energy-momentum 
conservation considerations is presented here. It is shown that the baseline of the two particle 
correlation function can be described using this technique. 



1. Introduction 

Since November 2009 the Large Hadron Collider (LHC) has produced 14M proton-proton 
collisions at ^=900 GeV in ALICE (A Large Ion Collider Experiment) [lJ.The ALICE detector 
is designed to investigate matter at the highest energy densities, where it is supposed to be 
deconflned. This state of matter is composed of free quarks and gluons, therefore known as the 
Quark Gluon Plasma (QGP), and it is believed that right after the Big Bang matter was in this 
state. The QCD phase diagram (Fig. [I]) shows the different phases of nuclear matter, and it is of 
great interest to study the properties of the QGP like the temperature of the phase transition, 
the critical point, the energy density, time evolution, lifetime of the system and flow-like effects. 
The present study is concerned with the measurement of the size of the particle emission region 
in proton proton collisions at ^=900 GeV using femtoscopy, and how it is affected by energy 
momentum conservation effects. 

2. Femtoscopy in ALICE 

Femtoscopy [2], or HBT [3J, is a technique that measures the shape and size of the particle 
emission region by measuring pairs of identical particles simultaneously, see Fig. [2j Pions are 
most commonly used in this analysis, and because they are indistinguishable (bosons) the two- 
particle wave function has to be symmetrized as follows: 
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Figure 1. The QCD phase 
diagram shows the different phases 
of nuclear matter. 




Particle 
Source 



P, 



'"""•''»».„ ****** Q=p -p 



Detector 



Figure 2. Two identical parti- 
cles are emitted and detected simul- 
taneously. There are two possible 
ways in which the particles could 
have been emitted. 



The probability of detecting two such particles is : 

|* 2 | - l + cos((pi -p 2 ) ■ (ri -r 2 )) = l + cos(Q- Ar), (2) 

where Q = p\ — p 2 is the momentum difference and Ar — T\ — r 2 is the relative emission point 
of the particles (Fig. [2|. Pairs with small relative momentum Q will be observed with higher 
rates because Bose-Einstein particles are likely to be found in the same state. 

For a source of particle emission p{r) a correlation function can be obtained by integrating 
over all possible pairs of pions resulting in the Fourier transform of the source function: 

C th (Q) = J P(Q)p(r 1 )p(r 2 )dr 1 dr 2 = 1 + p(Q) (3) 

In particular for a Gaussian source p(r) ~ exp _r * 2 / 2jR2 one obtains: 

C th (Q) = l + \e- R2 Q\ (4) 

where A is the correlation strength. This equation is used to fit the experimental correlation 
function, defined as the ratio of the two particle and single particle distributions: 
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The two particle distribution is obtained by correlating pions from the same event, and the single 
particle distribution is obtained from mixed events. The data shown in Fig. [3] is from 250 k 
proton-proton events at = 900 GeV recorded by ALICE in December 2009. The higher counts 
at low Q are the Bose-Einstein enhancement, but there are also long range correlations present, 
which are caused by non-femtoscopic effects such as jets and energy momentum conservation 
induced correlations (EMCICs for short). 



In Fig.Q C exp (Q) is plotted together with a Gaussian fit. It is evident that the simple Gaussian 
fit neither describes the low Q enhancement, nor the long range correlations. This implies that 
the particle emission source for proton-proton collisions is non Gaussian and other functional 
forms such as exponential or Lorentzian can be used. In this paper we keep using a Gaussian 
as we are interested in describing the baseline. 
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Figure 3. Two particle distribu- 
tions for real and mixed events. 



Figure 4. Correlation function 
with a Gaussian fit. 



3. Energy- momentum conservation induced correlations 

In order to determine the source size better it would be very useful to find a functional form for 
the sloped baseline observed in Fig.Q. A linear or quadratic slope can be used, but these are 
ad-hoc parameterizations, and it would be better if the baseline could be determined from first 
principles. It has been shown [I] that particle distributions are affected by EMCICs and if there 
are no other correlations present the two particle correlation function becomes: 
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where N is the event multiplicity. Equation [6] characterizes the EMCICs and it would be possible 
to remove their effect from the experimental correlation function. However, in experiment 
iV, (p^), (p^), (E) and (E 2 ) cannot be measured, simply because not every single particle is 
detected. These unknown quantities can be parametrized instead, and the resulting equation 
can be used as a fit to the data: 
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where the M parameters are defined as: Mi 
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N((E 2 )-(E) 2 ) anc ^ n °ta-tion {X} represents histograms of two-particle quantities that 
can be measured in experiment (see Fig [5]). 

Monte Carlo data, which does not have femtoscopic correlations but does guaranty energy 
and momentum conservation, has been used to test Eq. ([7| as a fitting function for the baseline. 
Figure [6] shows a good agreement between data and the fitted function. 

When all correlations are present the total correlation function can be written as: 



C(Q) = <S> fem to{Q) x C E mcic(Q) 



(8) 



where §f e mto = 1 + \e~ R2 ® 2 is the pure femtoscopic effect. Equation ^ is then used as the 
fitting equation to the real data. The results are shown in Fig. [7] and it can be seen that the 
baseline is described quite well by this method. The oscillations in the fitting function at large 
Q arise from the limited statistics. 



EMCIC histograms 



- <E1 +E2> 
<E1*E2> 
<Pt1*Pt2> 

- <Pz1*Pz2> 

- EMCICs 



ALICE Work in Progress 

ALICE 



^0 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Q(GeV/c) 



Figure 5. Th e different EMCIC 
histograms {X} are shown together 
with Cemcic{Q)- 





Figure 6. Eq. Q is used to fit 
Monte Carlo Data. 



Figure 7. Fit to the data using 
Eq. ®. 



4. Summary 

Two particle correlations for proton-proton collisions at 900 GeV exhibit long range correlations 
that make the determination of the baseline difficult. These correlations are primarily from 
energy momentum conservation and it has been shown here that they can be quantified and 
used to determine the baseline correctly. The EMCICs introduce more parameters into the 
fitting function, but the parameterization is obtained from first principles and these preliminary 
results look promising. Future work will be to separate the correlation function and EMCICs 
in pair average momentum bins (kr — Pl ~^ P2 ^j and fit all bins simultaneously with one set of 
parameters Mi through M4, thereby constraining more the parameter space. Three dimensional 
correlation functions in Q ou t, Qsidei Qlong an d its decomposition in spherical harmonics will also 
be studied. 
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